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Abstract 

We calculate the interaction potential between N atoms and NH molecules and 
use it to investigate cold and ultracold collisions important for sympathetic cool- 
ing. The ratio of elastic to inelastic cross sections is large over a wide range of 
collision energy and magnetic field for most isotopic combinations, so that sym- 
pathetic cooling of NH molecules by N atoms is a good prospect. However, there 
are important effects due to a p-wave shape resonance that may inhibit cooling in 
some cases. We show that scaling the reduced mass used in the collision is ap- 
proximately equivalent to scaling the interaction potential. We then explore the 
dependence of the scattering properties on the reduced mass and explain the reso- 
nant effects observed using angular-momentum-insensitive quantum defect theory. 

1 Introduction 

At temperatures below about 1 mK, atoms and molecules enter a fully quantal regime 
where their de Broglie wavelength is large compared to molecular dimensions. In 
this regime, collision cross sections and reaction rates are dominated by long-range 
forces 1,2 and resonance phenomena.— It is likely to be possible to control reaction 
rates by tuning scattering resonances with applied electric and magnetic he\ds.^~— At 
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even lower temperatures, below about 1 /iK, trapped atoms and molecules form quan- 
tum gases such as Bose-Einstein condensates and Fermi-degenerate gases, in which 
every molecule occupies the lowest allowed translational state in the trap. The quan- 
tum gas regime offers additional possibilities for a new form of quantum control, in 
which chemical transformations are carried out coherently on entire samples of ultra- 
cold atoms and molecules. 

There have been enormous advances towards these goals in the last few years. In 
particular, it is now possible to produce ultracold alkali metal dimers in their rovibronic 
ground states from ultracold atoms both by photoassociation 7 " 9 and by magnetoasso- 
ciation followed by stimulated Raman adiabatic passage (STIRAP). 1Q ' n However, for 
an alkali metal dimer even the ground rovibronic state has nuclear spin hyperfine struc- 
ture, 12 ' 13 and the resulting splittings are comparable to the kinetic energies involved in 
ultracold collisions. For the case of 40 K 87 Rb, microwave transitions have been used 
to transfer the ground-state molecules selectively between different hyperfine and Zee- 
man levels. 1Q ' 14 In a very recent development, Ospelkaus et ai— have studied reactive 
collisions of such state-selected molecules, both with one another and with ultracold 
K and Rb atoms. They observed remarkable selectivity of the resulting reactions, in 
which flipping the spin of a single nucleus could cause dramatic changes in the outcome 
of a collision.— 

Methods that form molecules from ultracold atoms can be applied only in cases 
where the atoms themselves can be cooled. In practice this restricts them to the alkali 
metals, the alkaline earths, and a few other elements. These species have a fairly lim- 
ited chemistry. In order to cool a wider class of molecules, including polyatomic ones, 
a number of direct cooling methods 17-21 have been established over the last decade. 
Among these methods, buffer-gas cooling is based on the particularly simple idea of 
cooling molecules by elastic collisions with cold He gas. If the molecules are param- 
agnetic and in low-field-seeking states, they can be confined in a magnetic trap. The 
temperatures which can be achieved in buffer gas cooling method are limited by the 
vapour pressure of the buffer gas (ca. 400 mK for 3 He), but the method is particu- 
larly valuable for two reasons: (i) it is very general and can in principle be applied 
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to any paramagnetic species, provided that a detection scheme is available; (ii) it can 
produce large numbers and high densities of cold molecules. Buffer-gas cooling has 
been reported for a variety of molecules including CaH, — CaF, 22 N H 23 ' 24 , CrH and 
MnH 25 and also for a number of paramagnetic atoms. 26-28 Buffer-gas cooling followed 
by evaporative cooling has recently been used to achieve Bose-Einstein condensation 
with no laser cooling for metastable helium. 29 

The direct methods established so far are limited to temperatures of 10 to 100 mK 
and above. To cool the molecules further, to the jiK regime, second-stage cooling 
methods must be developed. The most promising and conceptually the simplest method 
is sympathetic cooling, in which the molecules are cooled by collisions with an atomic 
gas that can itself be cooled to the ultracold regime, such as an alkali metal. The 
most robust trapping methods for molecules work for low-field-seeking states, which 
are never the lowest possible state in an applied field. Inelastic collisions can there- 
fore occur, and either heat the trapped system or eject the molecules from the trap. 
Sympathetic cooling can thus be successful only if elastic collisions dominate inelastic 
ones, and it is usually stated that the ratio of elastic to inelastic cross sections must 
be 100 or more. Sympathetic cooling was initially developed as a cooling method 
for trapped ions. 30 More recently it has been used to achieve sub-Kelvin temperatures 
for polyatomic ions 3 1 and has also been used to produce ultracold neutral atoms with 
scattering properties that are not suitable for evaporative cooling, such as 41 K. 32 

Sympathetic cooling for molecules has not yet been achieved, but several propos- 
als have been explored. It was initially proposed for NH molecules colliding with Rb 
atoms 33 and studied in more details by Lara et al. 34,35 for OH colliding with Rb. Both 
OH and NH molecules interact very strongly with Rb and the anisotropy of the in- 
teraction potential is large compared to the molecular rotational constant. The large 
anisotropy implies large couplings between channels with different n (monomer ro- 
tation angular momentum) and L (end-over-end angular momentum) quantum num- 
bers, and Lara et al. showed that this resulted in large inelastic cross-sections in the 
ultracold regime. The remedies they suggested to improve sympathetic cooling and 
decrease inelastic cross sections were: (i) to use light atoms as coolants, in order to 
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increase the heights of centrifugal barriers and suppress inelastic channels; (ii) to find 
atom-molecule system with much smaller anisotropy in the interaction potential. 

Soldan et al— considered the possibility of reducing the anisotropy by using alkaline- 
earth atoms (Ae) as collision partners for NH molecules. They showed that the neutral 
states of Ae-NH systems are coupled to ion-pair states Ae + NH , with crossings be- 
tween the neutral and ion-pair surfaces at linear geometries. For Sr and Ca atoms the 
crossings occurs at energies below the atom-molecule threshold, so will be accessible 
in low-energy collisions. However, for Be-NH and Mg-NH the crossings occurs at 
energies more than 1000 cm -1 above the atom-molecule threshold. In these systems, 
the ion-pair state is likely to be inaccessible, so it is reasonable to carry out collision 
calculations on a single covalent surface. In addition, the potential energy surface for 
Mg-NH turned out to be only weakly anisotropic. Wallis and Hutson 37 carried out 
quantum scattering calculations of spin relaxation collisions (in magnetic fields) and 
showed that sympathetic cooling of NH by collisions with Mg atoms should be achiev- 
able if the molecules can be precooled to about 10 mK. 

Sympathetic cooling has also been considered for NH3 and ND3. In this case the 
molecules are initially slowed in a Stark decelerator. 38 Zuchowski et al.— surveyed the 
interaction potentials for NH3 interacting with alkali-metal and alkaline-earth atoms. 
Zuchowski and Hutson 40 then carried out quantum scattering calculations on collisions 
of ND3 with Rb atoms and showed that molecules that are initially in the upper compo- 
nent of the ammonia inversion doublet are likely to undergo fast collisional relaxation 
to the ground state, and that this is likely to prevent sympathetic cooling of molecules 
trapped in low-field-seeking states in an electrostatic trap. 41 However, there is a good 
prospect for sympathetic cooling of ammonia molecules in high-field-seeking states, 
even with magnetically trapped atoms, because the terms in the hamiltonian that might 
cause spin-changing collisions of the Rb atoms are very small. High-field-seeking 
states of ND3 can be confined in an alternating current trap. 42 

Recently, Hummon et al.— demonstrated buffer-gas cooling and trapping of N ( 4 S) 
atoms and simultaneous co-trapping of NH molecules. Subsequent work 44 has demon- 
strated N atom densities around 5 x 10 12 cm~ 3 and lifetimes around 10 s. This offers 
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the possibility of cooling the atoms further with atomic evaporative cooling, which has 
already been achieved for metastable helium and Cr atoms. 45 i 46 

A gas of N atoms is potentially an excellent coolant for a sympathetic cooling 
experiment. The N atom has a very small polarizability compared to Group I and 
Group II elements and this results in low Cg coefficients and small anisotropies of 
the interaction potentials with molecules. The N atom also has a relatively low mass, 
which results in higher centrifugal barriers and stronger suppression of inelasticity for 
particles scattered with L > 0. 

This paper presents theoretical studies of cold and ultracold collisions of N atoms 
with NH molecules, in order to investigate the possibility of sympathetic cooling of NH 
by atomic nitrogen. Since the cross sections depend strongly on the reduced mass of the 
collision system, we consider four isotopic combinations of N-NH systems, with each 
of the two N atoms being either 14 N or 15 N. We assume that both N and NH are initially 
in their magnetically trappable spin-stretched states, with the maximum possible values 
of the electron spin projection numbers. For such states only spin relaxation (and not 
spin exchange) can occur and only the sextet interaction potential contributes. We 
report calculations of the sextet potential for N-NH and explore the behaviour of cross 
sections as a function of collision energy and magnetic field. We discuss the sensitivity 
of the scattering results with respect to uncertainties in the interaction potential. Finally, 
we analyze the behaviour of the shape resonances in terms of angular-momentum- 
insensitive quantum-defect theory (AQDT). 47 

2 N-NH potential 

The total spin of the N( 4 S) + NH( 3 £~) system can be |, | or j. The chemical reaction 
N + NH — >• N2 + H, which occurs principally on the doublet surface, has been studied in 
detail by Varandas and coworkers^"— and by Francombe and Nyman.— It was shown 
that the doublet N-NH system forms an N2H complex without a potential barrier along 
the minimum energy path. A very small barrier exists between the N2H complex and 
N2 + H products and overall the reaction of forming N2+H yields 6.33 eV of energy. 
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Table 1: Basis-set dependence of the N-NH interaction energy at the global minimum 
for F12 calculations. The complete basis set (CBS) extrapolation was obtained with the 
correlation energy functional E(X) = A +BX~ 357 where X is the maximum angular 
momentum of electronic basis set. 



basis set Ei nt (cm ) 



aug-cc-pVTZ 


-79.18 


aug-cc-pVQZ 


-86.03 


aug-cc-pV5Z 


-88.51 


aug-cc-pV6Z 


-89.33 


CBS 


-90.47 


F12 /aug-cc-pVTZ 


-89.10 



To our knowledge, no studies of quartet or sextet states of N-NH have been published, 
though we are aware of work in progress by Tscherbul and coworkers.— 

To obtain the sextet interaction potential we applied the recently developed ex- 
plicitly correlated, unrestricted coupled-cluster method with single, double and non- 
iterative triple excitations [UCCSD(T)].~~- 5 We used the aug-cc-pVTZ basis set of 
Peterson et al., 56 which is designed specifically for use with explicitly correlated cal- 
culations. The results from the explicitly correlated (F12) calculation are compared 
with those from UCCSD(T) calculations with uncorrelated basis sets in TableQ] it may 
be seen that the explicitly correlated approach dramatically reduces the error caused by 
using unsaturated basis sets. A fixed NH bond length of 1.0367 A was used in all the 
calculations. 

The potential energy surface was obtained by carrying out explicitly correlated 
UCCSD(T) calculations on a grid in Jacobi coordinates (Ri,9j), where R is the in- 
termolecular distance measured to the NH center of mass and 9 is the angle between 
the NH bond vector and the vector from the NH center of mass to the N atom. The ra- 
dial grid R{ was from 2.5 to 10 A in 0.25 A steps and the angular grid Oj was a set of 1 1 
Gauss-Lobatto quadrature points, which include the two linear geometries. All inter- 
action energies were corrected for basis-set superposition error using the counterpoise 
procedure. 58 

Radial interpolation is carried out using the reproducing kernel Hilbert space (RKHS) 
metho d 59 ! 60 to evaluate V (R, 0j) for arbitrary R and given 9j. At each distance R, the 
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Table 2: Van der Waals coefficients for N-NH (Eha'o) from density-functional calcula- 
tions. 

n, A C„x 

6.0 33.50 

6.2 11.44 

7.1 60.87 

7.3 55.28 
8,0 717.14 

8.2 988.87 



potential is expanded in Legendre polynomials Px (cos 9) for A up to 8, 



V(/?,0)=£V A (/?)f A (cos0). (1) 
X 

The coefficients Vx (R) are obtained by integrating the ab initio potential using Gauss- 
Lobatto quadrature. 35 

To provide an improved description of the long-range interaction, we impose an 
analytical representation on the long-range part of the components of the projected 
potential, 

V%(R) = -t ic n ,xR- n - (2) 

n=6X=0 

The Van der Waals coefficients are given in Table [2] and were calculated with the 
restricted open-shell coupled Kohn-Sham method 61 with asymptotically corrected 62 
PBEO functional. 63 We connect the long-range function smoothly to the supermolecu- 
lar potential using the switching function— 



/w4+r in ?( 3 - sin2 ?)' (3) 



where x = R ~ b h t*~ a with a = 7 A and b = 1 1 A. f (R) = for R < 7 A and f(R) = 1 
for R > 1 1 A. 

The potential energy surface for N-NH is shown in Fig. Q] It has two minima of 
comparable depths at linear geometries: 89.1 cm -1 at N-NH and 76.4 cm -1 at N- 
HN. The two minima are separated by a saddle point near the T-shaped geometry. The 
anisotropy of the potential near the Van der Waals minimum is about 40 cm -1 , and the 
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Table 3: Characteristic points on the N-NH potential energy surface. Energies are 
given in cm -1 , R in A. 



Global minimum 


Saddle point 


Secondary minimum 


R,6 3.70,0 


3.76, 92° 


3.49, 180° 


-89.1 


-39.2 


-76.4 
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Figure 1 : The ab initio interaction potential of N-NH. Contours are labeled in cm" 1 . 
The angle 0=0 corresponds to the N-N-H geometry. 

dominant contribution to the anisotropy arises from V2(R). 
2.0.1 Interaction potential uncertainty 

An important problem in electronic structure theory is the estimation of error bounds 
for calculated interaction energies. Since scattering calculations at very low energies 
depend strongly on the details of the interaction, in this section we discuss the uncer- 
tainty of the N-NH interaction potential obtained here. 

The largest contributions to the uncertainty of the interaction potential arise from 
the approximate treatment of electronic correlation and the incompleteness of the elec- 
tronic basis set. We expect that the effect of the neglecting vibrations of the NH 
molecule is much less important, as are relativistic and nonadiabatic effects. 

First we need to estimate how well the UCCSD(T) method works for sextet N-NH. 
To explore this, we performed 7-electron full configuration-interaction (FCI) calcula- 
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tions of the interaction energy. We included all electrons arising from the H atom and 
the 2p electrons of N atoms. With FCI it was possible to use only a very small basis set 
(6-3 1G, augmented with spd midbond functions with an exponent 0.4). For this small 
basis set, we compared the contribution to the correlation part of the interaction energy 
(i.e., the intramonomer correlation and the dispersion energy) with the UCCSD(T) re- 
sults for several linear geometries N-NH and N-HN . The FCI correlation energy to be 
larger than the UCCSD(T) correlation energy by 1 to 1 .5%, for a wide range of R and at 
both linear geometries. To a good approximation we expect the ratio E^. / Ec^ SD ^ 
to be constant in different basis sets. This suggests that the global minimum energy 
obtained with the coupled-cluster method is underestimated by ca. 1.5 cm -1 . 

The basis set convergence pattern shown in Table Q] yields a complete basis-set 
limit of the global minimum depth of 90.47 cm - '. This is 1.37 cm" 1 more than in 
the method used for the complete surface here. We also performed test calculations 
including additional core-valence basis functions that are absent in the basis set used for 
the complete surface potential. The interaction energy at the global minimum obtained 
with aug-cc-pCVTZ is approximately 1 cm smaller than for basis sets with no core- 
valence functions. 

In summary we can set the error bounds on the interaction potential at the global 
minimum between —1 and +3 cm" 1 , which is approximately between —1% and +3%. 

3 N-NH scattering calculations 

The Hamiltonian of the NH molecule may be written 

#nh = b NH N 2 + 7 N-S 

A ss £(-iry 2 ,_ ? (f)[S®5f). (4) 
9 

The three terms are, respectively, the rigid rotor Hamiltonian, the spin-rotation interac- 
tion and the intramonomer spin-spin interaction. The numerical values of the constants 
used in the present work are Z?nh = 16.343 cm" 1 ,— y = —0.055 cm" 1 and Ass = 0.92 
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cm 1 . 66 The NH molecule is assumed to be in its ground vibrational state. 

The Hamiltonian of the N-NH collision system in a magnetic field may be written 

H = - Sr7& r + 2^ +//NH +//z + Vss +VM 0) - (5) 

Here L? is the operator for the end-over-end angular momentum of N and NH about 
one another, Hz represents the Zeeman interaction of N and NH with the magnetic 
field, Vss is me (anisotropic) intermolecular spin-spin interaction, and Vi nt (R, 0) is the 
intermolecular potential. 

The convention for quantum numbers in this paper is as follows: L and Mi denote 
the end-over-end angular momentum and its projection onto the space-fixed Z axis 
defined by the magnetic field. Monomer quantum numbers are indicated with lower- 
case letters to avoid confusion with those of the collision system as a whole. The 
spins and spin projections of the N and NH molecules are denoted by sq and m S A, 
and m S Q, respectively. The rotational quantum number of the NH molecule and its 
projection are denoted uq and m n g The projection of the total angular momentum, 

M tot =M L + m nB + m S B + m sA , (6) 

is rigorously conserved in a collision, but the total angular momentum itself is not, 
except at zero field. It is convenient to carry out scattering calculations is a fully un- 
coupled basis set, |sawi s a) I-sbtWjb) \ n a m na) \J-Ml) ■ We have written a plug-in routine for 
the MOLSCAT scattering program, 67 implementing all the matrix elements required 
for scattering calculations in this basis set. 
set. 

The total spin S of a system made up of an open-shell atom and an open-shell 
molecule can take values between |% ~ $b I an d $A + s b ■ For N-NH the allowed values 
are S = j, | and |, corresponding to doublet, quartet and sextet, respectively. The 
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interaction potential Vi nt (R, 9) may be written in terms of projection operators, 
V iat (R,8)= jf \S)V S (R,9)(S\ 

S=-|-«A+-«B 

and the general matrix element of Vi n t(^, 0) in our basis set is 



The three interaction potentials Vs(R,9) differ only by short-range Pauli exchange 
terms. They have the same long-range coefficients, so become degenerate once the 
N atom and NH molecule are far enough apart that their valence shells do not over- 
lap. The doublet surface has a potential well several hundred times deeper than the 
Van der Waals sextet state, so that full quantum calculations including the doublet 
state would require very large basis sets of rotational functions and could not be con- 
verged. In the present work we therefore approximate the operator V mt (R, 0) operator 
by taking V$ = V5/2 for all spin states. This approximation is legitimate because we 
are primarily interested in N-NH collisions between magnetically trapped atoms and 
molecules, with m s x = s a = § and m S Q = sq = 1. These are spin-stretched states, 
and V 3 n and V\ n have no matrix elements (diagonal or off-diagonal) involving spin- 
stretched states. When this approximation is made, orthogonality relations for the 3 j 
symbols reduce Eq. [8] to a form diagonal both in m sA and m s g. The explicit expres- 
sion for (rigm n gLMi\Vs(R, 9)\ngm' nB L' M' L ) is the same as for scattering of NH from a 
closed-shell atom, 5 with the addition of factors 8 m m i . 



(s A m sA SBm S BnBm nB LM L \V inl (R, 0)|s A ^ A s B 'niB«B' M «B L ' M L> = 
) 2.s A+2iB -m lA -m lB -M i ^ + ^ ^ miiRLML \y s ^ d)\n B m' nB L> 'M' L ) 




(8) 
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The intermolecular spin-spin interaction has matrix elements 



(s A m sA s B m sB n B m nB LM L \Vss\ s A m' sA s B m sB ri B m nB L ' M' L ) 



V36X(R)$n B n>J m , lB ,„> B (~lY A+s »- m * A - m *- ML lsA(sA + l)(2s A + l)s B (s B + 1)(2, B + 1)(2L+ 1)(2L' + 



The spin-spin coupling constant A(/?) is E^a 2 a^/R? , where is the Hartree energy 
and a is the fine-structure constant. 

The matrix elements for NH monomer operators are the same as for scattering of 
NH from a closed-shell atom, 5 with the addition of factors 8,„ ,„/ . 

If one or both of the colliding species is not in a spin-stretched state (with the high- 
est possible value of ms), the system will undergo very fast spin exchange driven by the 
difference between the S = |, | and S = j potentials. For spin-stretched states, how- 
ever, spin exchange cannot occur and only spin relaxation is possible. There are two 
mechanisms for spin relaxation. The first is similar to the well-known mechanism of 
spin relaxation for spin-stretched states of alkali metal atoms, and arises through direct 
coupling of the initial state m sA = +j,m sB = +1 (with n B = m„ B = 0), to final states 
with m sA and/or m sB reduced by 1 by the intermolecular spin-spin interaction term and 
Ml increased to conserve M m . Such transitions are relatively slow, because the in- 
termolecular spin-spin interaction is weak. The second mechanism is that described 
by Krems and Dalgarno. 68 For the same initial state, the intramonomer spin-spin in- 
teraction mixes n B = with n B = 2, and even in a magnetic field it mixes m„ B = 
with m nB = ±1,±2. The states with m nB = +1.+2 have m sB =0,-1. The states with 
m„ = 1 , 2 are then coupled by the anisotropy of the interaction potential to n = 0, m n = 
states with changed Mi but the same m B (which is lower than in the initial state). This 
mechanism is also expected to be fairly weak for a low-anisotropy system such as N- 
NH: the n = and n — 2 rotational levels of NH differ in energy by 96 cm" 1 , while the 
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potential anisotropy V*2(i?) that couples them is a short-range interaction that is never 
greater than 40 cm -1 in the energetically accessible region. For both mechanisms, 
spin relaxation is suppressed for s-wave scattering (L = 0,Mi = 0) at low energies and 
fields because the conservation of M tot requires M' L ^ and therefore L' > 0, producing 
a centrifugal barrier in the outgoing channel. 

We carry out scattering calculations with the MOLSCAT package. 67 The coupled 
equations are solved using the hybrid log-derivative/ Airy propagator of Alexander and 
Manolopoulos. 69 We used the fixed-step log-derivative propagator from 2.8 A to 70 
A with an interval size of 0.08 A, followed by a variable-step Airy propagation out 
to 400 A. We carried out convergence tests on state-to-state cross sections both in the 
s-wave regime and at energies up to E = 1 K, at fields of B = 200 G, 1000 G and 
2 T. In all cases a basis set with n = . . . 3 and L = . . .7 gave convergence to within 
approximately 1% for all state-to-state cross sections. This basis set was therefore used 
in all the remaining calculations. 

4 Results 

Fig. |2] shows the Zeeman energy levels of the noninteracting N+NH system. In the 
buffer-gas cooling experiment, 43 both atoms and molecules are trapped in their low- 
field-seeking state with m S A = \ and «b = w«b = 0, m S B = 1 ■ The experiment has 
already achieved temperatures around 550 mK, and at this temperature atoms with an 
energy of 5kT sample magnetic fields up to 2 T in a quadrupole trap. However, as the 
temperature decreases, so too will the magnetic fields sampled. We therefore consider 
collision energies from 10 jiK to 1 K and fields from 10 G to 2 T 

The N atom is considerably less polarizable than alkali metal or alkaline earth 
atoms. As a result, the dispersion coefficient C^ q for N-NH is considerably lower 
than for most metal atom - molecule systems that have been considered previously as 
candidates for sympathetic cooling. Together with a low reduced mass, this results in 
relatively high centrifugal barriers for L > partial waves: 14 mK for L = 1,71 mK for 
L = 2, 120 mK for L = 3, etc. The high centrifugal barriers also mean that quite small 
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m sA =3/2m sB =l 

m sA =3/2 m sB =0,m sA =l/2 m sB =l 

m sA =l/2 m sB =0,m sA =-l/2 m sB =l 
m sA =3/2m sB =-l 

m sA =-l/2 m sB =0,m sA =-3/2 m sB =l 
m sA = l/2m sB =-l 

m sA =-3/2 m sB =0,m sA =-l/2 m sB =-l 
m sA =-3/2m sB =-l 



B I Gauss 

Figure 2: Energy levels of the noninteracting N+NH system in a magnetic field. The 
dotted red lines show the energy obtained by adding the d-wave centrifugal barrier 
height (71 mK) to the levels with M tot = | and M tot = |. The crossings between the 
red lines and the initial-state energy indicate the fields above which s-wave inelastic 
cross sections are no longer suppressed by centrifugal barriers. 

number of partial waves are needed to converge cross sections: for example, including 
contributions from L up to 4 is sufficient to obtain convergence up to about 0.5 K. 

4.1 Close-coupling calculations 

Calculated elastic and inelastic cross sections for the four different isotopic combina- 
tions are shown as a function of collision energy E in Fig. [3] for representative magnetic 
fields of 200 G, 300 G, 1000 G and 2 T. In a simple hard-sphere model of sympathetic 
cooling, neglecting inelastic collisions, the temperature relaxes towards equilibrium 
and reaches a lie point after [m\ + m-i) 1 / , 2mini2 collisions, 70 where m\ and mq. are the 
masses of the two species. For sympathetic cooling to be successful we need the ratio 
of elastic to inelastic cross sections to be much larger than this. The calculated ratios 
are shown in Fig. [4] for the most part they are more than 50 at collision energies above 
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Figure 3: Elastic and inelastic cross sections for N-NH scattering for N- NH, N- 
15 NH, 15 N- 14 NH and 15 N- 15 NH for different magnetic fields. The elastic cross section 
(black line) is almost independent of field. The positions of the p- and d-wave exit- 
channel barriers are marked with vertical lines. 
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Figure 4: The ratio of elastic to inelastic cross sections for 14 N - 14 NH, 14 N - 15 NH, 
15 N - 14 NH and 15 N - 15 NH systems for different magnetic fields. The value of 100 
typically required for sympathetic cooling is marked with a dotted horizontal line. 



about 1 mK, indicating that sympathetic cooling of NH by N is likely to be feasible. 

Several different effects are evident in Fig. [3] The first is that the cross sections 
enter the s-wave regime, where they are proportional to £ -1 / 2 , at quite different en- 
ergies for different isotopic species. This occurs because of p-wave resonant effects. 
Once in the s-wave regime, however, the inelastic cross sections generally decrease at 
magnetic fields below about 300 G, because of the centrifugal barriers in the outgoing 
channels. Since atoms and molecules in a quadrupole trap sample lower and lower 
fields as the temperature is decreased, this indicates that sympathetic cooling will be- 
come increasingly effective as the temperature is lowered, as predicted for Mg-NH. 37 
Lastly, inelastic collisions are also suppressed for very high magnetic fields. All these 
effects will be discussed in more detail below. 

Spin relaxation collisions can change m S A for the N atom, m s g for the NH molecule, 
or both. Fig.|5]shows the state-to-state cross sections for the most important final states 
for 15 N- 15 NH as a function of energy at two different fields. It may be seen that domi- 
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Figure 5: State-to-state inelastic cross sections for N - NH for weak (10 G) and 
strong (0.1 T) magnetic fields. 
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nant final states are those in which m S A and/or m S Q has changed by 1 . These collisions 
are driven by the intermolecular spin-spin interaction. Transitions that change m S Q by 
2 can occur only by the second mechanism described in Section[3]above, involving the 
potential anisotropy, and are seen to be very much weaker except in a small resonant 
region. 

The 14 N- 14 NH system shows behaviour quite different from the others, with a large 
peak in the inelastic cross sections near 10 mK which reduces the elastic-to-inelastic 
ratio to around 10. This ratio may not be high enough for effective sympathetic cooling 
from an initial temperature of tens of milliKelvin. The peak appears at the same energy 
for all values of the field. It arises from a p-wave shape resonance in the incoming 
channel, as discussed in section |431 below. For the larger reduced masses of the other 
isotopic combinations, the quasibound state responsible for the shape resonance drops 
below threshold and becomes a true bound state. Thus the other isotopic combinations 
do not exhibit this feature and have more favourable properties for sympathetic cooling. 

The 15 N- 15 NH system exhibits d-wave shape resonances for collision energies of 
50 to 70 mK, but they are much weaker than the p-wave resonance for 14 N- 14 NH and 
their presence does not strongly affect the total inelastic cross section. 

The L = 2 centrifugal barrier plays a crucial role in spin relaxation in the ultracold 
regime. For an incoming channel with L = 0, spin relaxation requires outgoing L > 
2. If the energy difference between the incoming and outgoing channels is smaller 
than the height of the L = 2 centrifugal barrier, the s-wave inelastic cross section is 
suppressed (see Fig.[2j». The s-wave state-to-state cross sections are shown as a function 
of magnetic field in Fig.|6]for 14 N- 14 NH at a collision energy of 50 jJ.K. The inelastic 
cross sections generally decrease at magnetic fields below about 500 G, though there 
is a dip in each state-to-state component between 100 and 300 G. These dips are due 
to suppression of the inelastic cross sections in the wings of resonances, as described 
by Hutson et al. ;— in this case the resonances concerned are shape resonances in the 
d-wave outgoing channels. 

The suppression of inelastic scattering by the centrifugal barrier is clear in the total 
inelastic cross section only for systems with reduced mass larger than for 14 N- 14 NH. 
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Figure 6: Inelastic s-wave cross section as a function of magnetic field for 14 N- 14 NH 
at a collision energy of 50 jiK. The dashed vertical lines correspond to the fields for 
which the kinetic energy release from M tot = | to M lot = | and M tot = A states are 
equal to the height of the d-wave barrier. 

For the 14 N- 14 NH system itself, the p-wave contribution is very strong even at very low 
energies. In fact, the p-wave enhancement of the inelastic cross section between 50 jiK 
and 1 mK is so strong that the total inelastic cross section does not follow the E~T- 
power-law dependence expected from the Wigner threshold laws at these energies. 

Suppression of inelastic collisions due to barriers in the outgoing channels de- 
creases as the magnetic field increases (so that the kinetic energy release increases). 
Eventually, however, the inelastic cross section reaches a maximum and starts to de- 
crease again. This occurs for all partial waves, and the total cross sections at a field of 
2 T are typically reduced by a factor of about 10 from their values at 0.1 T. Since for 
some isotopic combinations the low-field ratio of elastic to inelastic cross sections may 
not be large enough at temperatures of 1 to 10 mK, the application of a strong bias field 
to a magnetic trap offers a possible way way to improve the ratio. 

Suppression of inelastic cross sections at high fields has been observed for 0( 3 P)- 
He collisions, 71 for OH-OH 72 and for collisions of Cr atoms. 73 For small inelasticity, 
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Figure 7: The integrand of the distorted-wave Born approximation (Eq.fTOb for 14 N- 
14 NH at a collision energy of 1 mK for different kinetic energy releases. The incoming 
wavefunction is calculated for L = and the outgoing wavefunction for L = 2. 

the distorted-wave Born approximation gives 74 



where y/,- and \j/f are energy-normalized wavefunctions in the initial and final channels, 
Uif is the coupling between the channels, and kj is the wave vector in the incoming 
channel. Fig. [7] shows the integrand of Eq. [I0]for the intermolecular spin-spin term 
(R~ 3 ) at kinetic energy releases of 0.5 K, 1 K and 3 K, corresponding to fields of 
3800 G, 7600 G and 2.3 T for transitions with Am S A + Amjb = — 1 . It may be seen that 
there is significant oscillatory cancellation in the integral at high fields, when y/;(/?) 
and yff(R) oscillate out of phase with one another in the interaction region, and this 
combines with the effect of the resonances in the d-wave outgoing channels to produce 
the maxima in Fig. [6] The oscillatory cancellation occurs for arbitrary partial waves. 
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Figure 8: Elastic and total inelastic cross sections for 14 N- 14 NH at E = 5 mK and 
B = 50 G as a function of both reduced mass and potential scaling factors. 

4.2 Dependence on interaction potential and reduced mass 

Because of the crucial role played by the p-wave shape resonance for 14 N- 14 NH at 
energies up to 10 mK, it is important to investigate the influence of uncertainties in the 
interaction potential on the cross sections. The shape of the 2D interaction potential is 
complicated and the scattering properties might in principle depend on many param- 
eters. However, Gribakin and Flambaum 75 showed that for single-channel scattering 
the scattering length a behaves as 



a = a 1 - Ian ( <3> ^ 



(ID 



where, for a potential with long-range form —C(,R 6 , the mean scattering length a is 
0.956(2^C 6 //i 2 )5 and 



*= / (2iiV mt (R)/h 2 y dR 



(12) 



Although N-NH is a many-channel scattering problem, it is elastically dominated and 
Eq.Q~2]with V[ nt (R) replaced by Vq(R) reproduces the major features of the elastic scat- 
tering. Thus scaling fi is approximately equivalent to scaling the entire interaction 
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potential, and either of these scalings provides a good way to explore the variation of 
scattering length as a function of potential. Fig. [8] shows the elastic and total inelastic 
cross sections for 14 N- 14 NH at E = 5 mK and B = 50 G as a function of both reduced 
mass and potential scaling factors. It may be seen that the two scalings have a very 
similar effect, apart from a small shift in the Feshbach resonance around A = 1.08, 
which comes from a rotationally excited state of NH. 

We estimate the bounds on the accuracy of our potential to be between —1 and 
+3% of the well depth. We have therefore carried out calculations of the cross sections 
as a function of a variable reduced mass, parameterized by scaling factor jj. — >• Xfl for 
a collision energy of 5 mK. A scaling factor A = 1 corresponds to the reduced mass 
for 14 N- 14 NH. The result is shown in Fig. [9] There is a strong maximum in the total 
inelastic cross section near A = 1.012, due to a p-wave shape resonance in the elastic 
channel. A change of 1.2% in the potential is within the estimated error bound of our 
calculations. Enhancement of the cross sections due to the p-wave resonance might 
thus occur for heavier isotopic combinations than 14 N- 14 NH if our potential is slightly 
too deep. Although we believe it is more likely that our potential is too shallow than 
too deep, this cannot be ruled out. However, it is quite unlikely that enhancement due 
to the p-wave resonance would occur for the heaviest system, 15 N- 15 NH. 

4.3 AQDT analysis of shape resonances 

In this section we consider the N-NH scattering in the context of angular-momentum- 
insensitive quantum defect theory (AQDT).— The upper part of Fig. [10] shows the 
positions of (quasi)bound states for L = ... 2 as a function of the reduced-mass scaling 
factor A for values between 0.8 and 1.2. The L = bound state crosses the threshold 
at a value of the reduced mass much smaller than that for 14 N- 14 NH, well outside the 
estimated error bounds for the potential. There is thus no s-wave resonance in the 
scattering for any of the systems considered here. 

The p-wave shape resonance in the cross sections for 14 N- 14 NH arises from the 
quasibound state with L = 1, which is 6 mK above threshold for A = 1. As the reduced 
mass increases above this, the L = 1 bound state crosses the threshold at A = 1 .024 and 
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Figure 9: Elastic and inelastic cross sections for a collision energy of 5 mK as a function 
of scaled reduced mass. The dashed vertical lines indicate the values of jJ. for the 
systems concerned here, from left to right, 14 N- 14 NH, 14 N- 15 NH, 15 N- 14 NH and 15 N- 
15 NH, respectively. The solid vertical line shows the value A = 1.024 for which the 
scattering length a = 2d. 
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becomes a "real" bound state. This explains why we see no p-wave shape resonances 
in the cross sections for collisions with reduced mass larger than 1.024/Xi4.i4, i.e. for 
all systems containing at least one 15 N atom. For large values of fi the L = 2 bound 
state comes close to threshold, and this results in the (small) d-wave shape resonance 
that can be seen in the 15 N- 15 NH cross sections. 

The lower panel of Fig. [10] shows the s-wave scattering length as a function of 
scaling factor X. As expected, there is a pole near X — 0.9 as the L = bound state 
crosses the threshold. Within the estimated error bounds of the potential, a varies 
between 22.3 and 15.2 A, so that the elastic cross section for small E co \\ (in the Wigner 
regime) varies between 6250 and 2900 A 2 . 

In angular-momentum-insensitive quantum defect theory, 47 ■ 76 - 77 the scattering prop- 
erties of a system for arbitrary L can be predicted from only a few parameters: the 
s-wave scattering length a, the dispersion coefficient Ce and the reduced mass fi (and 
thus a). The positions where L > bound states cross threshold, and hence produce 
shape resonances, depend only on the relationship between a and a. In particular, when 
a = 2a there is an L = 1 bound state exactly at threshold, and systems with a slightly 
larger than 2a have a p-wave shape resonance at a collision energies below the height 
of the p-wave centrifugal barrier. This is the case for the 14 N- 14 NH system here. For 
a = a there is an L = 2 bound state exactly at threshold. The scattering length for 15 N- 
15 NH is 1.54a, which is close enough above a to produce a d-wave shape resonance at 
finite energy. The energies at which the p-wave and d-wave resonances appear can be 
read off the L = 1 and 2 lines in Fig.fTOl 

It should be noted that a change in the interaction potential would result in "slid- 
ing" the vertical lines representing the four isotopic combinations horizontally along 
Fig. [10] A range of behaviour can exist for potentials within our uncertainties of our 
calculations. Our calculations thus do not definitively identify which characteristics 
will be be observed for a particular isotopic combination. 

To study the shape of the resonant features more quantitatively, we can consider 
the Cl(E) functions introduced by Mies. 78,79 These functions give the connection be- 
tween a semiclassical JWKB description of scattering states (valid at large collision 
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Figure 10: Near-threshold bound states of sextet N-NH for L — ... 2 as a function of 
fl (upper panel), and the /x-dependence of the s-wave scattering length (lower panel). 
Field-free calculations with a structureless atom and diatom were used to obtain the 
bound states here. The dashed vertical lines indicate the values of fi for 14 N- 14 NH, 
14 N- 15 NH, 15 N- 14 NH and 15 N- 15 NH, respectively, from left to right. The red dotted 
line on the lower panel indicates 2d and the green dotted line indicates a. 
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Figure 1 1 : The p-wave transmission function Cj~ 1 (E) for N-NH systems with reduced 
masses corresponding to: 14 N- 14 NH (a/a = 2.34), 14 N- 15 NH (a/a = 1.94), and an 
artificial system with A = 1.031, for which a/a — > 2+. 

energies) and the near-threshold behaviour. The function C£ l (E) can be viewed as an 
enhancement factor in the short-range part of the wavefunction due to the presence of 
the long-range potential, including any resonant effects in the incoming channel. In 
Fig. QT|we show the C^ l [E) functions for p-wave scattering with A = 1, 1.023 and 
1.031 (the last of these values corresponding to the reduced mass of 14 N- 15 NH). As 
the scattering length decreases and reaches 2a, the height of the peak in C^ l (E) goes 
to +°°, as shown in Fig.[l2j and the energy at which the peak occurs approaches zero. 
The width of the resonance decreases, corresponding to increasing the lifetime of the 
quasibound state. The intensity of the resonance rapidly decreases once the scattering 
length is larger than 2a, corresponding to a bound state below threshold. 

5 Conclusions 

We have calculated a potential energy surface for N atoms ( 4 S) interacting with NH 
molecules ( 3 £~) in the spin-| (sextet) state, using unrestricted coupled-cluster calcula- 
tions with an explicitly correlated basis set. This is the surface that governs collisions 
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Figure 12: The height of the peak in the Cj~ 1 (E) function as a function of reduced-mass 
scaling factor X. 

of cold N atoms and NH molecules in a magnetic trap. We have used the surface to 
carry out quantum scattering calculations of cold collisions for different isotopic com- 
binations of N and NH, as a function of collision energy and magnetic field. 

The sextet potential energy surface is weakly anisotropic, with an anisotropy of ap- 
proximately 40 cm"' in the well region. The anisotropy is dominated by the P2(cos 8) 
Legendre, which mixes states with An = ±2 in the NH rotational quantum n. Since 
the anisotropy is smaller than the separation between the n = and 2 states, it causes 
relatively weak mixing during collisions and the scattering is generally elastically dom- 
inated. The inelastic cross sections are suppressed both at low energy and low field (by 
centrifugal barrier in the exit channels) and at very high field (by oscillatory cancella- 
tion due to the large kinetic energy release). For most isotopic combinations the ratio of 
elastic to inelastic cross sections is high enough, over a wide enough range of collision 
energy and magnetic field, that sympathetic cooling of NH by N is a good prospect. 

We have shown that a scaling of the interaction potential is approximately equiva- 
lent in its effects to a scaling of the collision reduced mass. We estimate our interaction 
potential to be accurate to within 3%. We have investigated scaling the reduced mass 
by up to 20% from the value for 14 N- 14 NH. The scaling revealed that there are ma- 
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jor effects arising for a p-wave shape resonance, which produce enhanced inelastic 
scattering for 14 N- 14 NH at low energies on our best potential surface. We have used 
angular-momentum-insensitive quantum defect theory (AQDT) to understand how the 
results change for different isotopic combinations. 

Scaling the potential energy surface, or equivalently the reduced mass, is a very 
useful tool for understanding cold collision calculations. In combination with AQDT, 
it can provide powerful insights into low-energy scattering for low-energy collisions 
where where only a few partial waves contribute to the scattering. 
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